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PREFACE 


A fortuitous set of circumstances brought together a 
group of people with diverse interests to address the 
important problem of finding ground water. It appeared 
that techniques developed for electromagnetic probing of 
the lunar interior, and techniques developed for the 
generation of high power audio frequencies, could be combined 
to make practical a magnetic inductive coupling system for 
the rapid measurement of ground conductivity profiles which 
are very helpful when prospecting for the presence and 
quality of subsurface water. 

Thus, we undertook the study and conceptual design of 
the system described herein which involves the measurement 
of the direction, intensity, and time phase of the magnetic 
field observed near the surface of the earth at a distance 
from a horizontal coil energized so as to create a field 
that penetrates the earth. From such observations, we 
planned to deduce the conductivity and stratification of 
the subsurface. 

As a result of our theoretical studies and a rudimentary 
experiment in an arid region we now can show that the approach 
is conceptually valid and that this geophysical prospecting 
technique deserves to be developed into a pragmatic system 
for the economical exploration of subterranean water 


resources . 
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1 . 0 Summary 

A geophysical probing technique, which we call SWEEP, 
based on a measurement of the inductive coupling between two 
coils deployed on, or near, the surface of the earth has 
been developed as an extension of earlier work. One coil 
is energized with an alternating current and the other is 
used as a means of measuring the apparent magnetic field at 
a distance. Because the system requires no physical contact 
with the ground, and because many different frequencies are 
used, the inductive technique avoids many of the problems, 
and resolves many of the ambiguities, inherent in other 
geophysical probing techniques. 

The SWEEP inductive coupling technique shows promise as 
a means whereby subsurface conductivity can be measured 
quickly, easily, and inexpensively. And, because subterranean 
water causes contrasts in subsurface conductivity, this 
technique shows promise as a useful tool for the detection 
of ground water in arid regions. 

Using a collection of available equipment, a field 
trial was conducted in an arid region near El Paso, Texas, 
to test the validity of the concept and the feasibility of 
making an automated system for the measurement of ground 
conductivity. In this field trial, a coil comprising twelve 
turns of #10 wire was laid on the ground in the form of a 
square with a side of 25 meters. This coil was energized 
with alternating current in the audio frequency range - with 
various pov.’er sources ranging from three kilowatts at the 
lowest frequency of 60 hertz to ten milliwatts at the 
highest frequency of 20 kilohertz. The resultant magnetic 
field was measured at distances from the source of 20 
meters to 500 meters, using a manually transported coil 
having 1200 turns with a cross-sectional area of one square 
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meter; this coil being arranged so that it could be rotated 
mechanically to either a horizontal or vertical position. 

A phase locked loop voltmeter was used to measure the 
magnitude and phase of the observed vertical and horizontal 
magnetic field components. 

Using extensions of prior theoretical analyses and a 
computer implemented model, we have been able to achieve an 
interpretation of the experimental data that agrees well 
with geological conditions of the field sites determined 
previously by much more laborious techniques. Specifically 
our interpretation yields an indication that at one test 
site there is a water bearing stratum nine meters below the 
surface of the earth, and at the other test site there is a 
water bearing stratum at a depth of 100 meters. These 
results agree with the depths at which water is found in 
nearby existing wells. 

To provide the feasibility demonstration described in 
this report, we used crude hardware, considerable manual 
labor during the field experiment, and lengthy 
calculations utilizing a large computer facility. However, 
with the knowledge gained through the development of our 
measurement and interpretation techniques, we feel that no 
further major technological advances are required to 
implement an automated system that would permit the 
economical mapping of subsurface conditions ir. arid regions 
at depths less than a few hundred meters. 

Conceptually, the ultimate SWEEP system would comprise 
two small instrumented vehicles, one carrying a horizontal 
source coil and a high power audio frequency generator, 
the other carrying a receiving coil and a coherent detector 
Data could be interpreted in real time in the field, and 
since there is no requirement for physical contact with the 
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ground, both vehicles could be maneuvered to permit rapid 
and economical mapping of large areas. The feasibility 
study and experiments reported here demonstrate that this 
concept is viable and realistic. 



2.0 Introduction 


The geophysical prospecting technique reported here 
involves the measurement of the complex inductive coupling 
between two coils deployed on or near the surface of the 
earch. A source coil, deployed horizontally, is energized 
with currents of the order of several amperes at frequencies 
ranging from a few hertz to a few kilohertz. A detecting 
coil, oriented alternately horizontally and vertically, at 
distances up to several hundred meters, is used in the 
measurement of the magnitude and phase of the field created 
just above the surface of the earth. Inversion of these 
data then yields a measure of the subsurface conductivity. 

Like other electrical prospecting methods, the SWEEP 
procedure measures electrical conductivity as a function of 
depth. And, because electrical conductivity of rocks and 
sediments depends strongly on the porosity, the amount of 
saturation, and the electrical conductivity of the 
pore-fluids, this conductivity data may be interpreted in 
terms of the presence or absence of water. 

The low frequency inductive technique promises several 
major advantages over other electrical prospecting techniques. 
First, a properly instrumented system is likely to be much 
faster to operate than existing techniques because there is 
no necessity for electrical contact with the ground and data 
can be collected while an instrumented vehicle is moving. 
Second, the alternating current frequency, which is one of 
the prime system variables, can be swept automatically so 
that a large amount of data can be collected rapidly without 
manual intervention. And, third, the problems associated 
with highly conductive layers near the surface, which may 
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be a severe handicap with other techniques, can be overcome 
by using a set of frequencies such that the surface layer is 
essentially transparent but the underlying strata can be 
delineated clearly. 

Further, and importantly, interpretation of data 
collected with the low frequency inductive technique 
promises to be less ambiguous than interpretation of data 
collected with existing techniques. In common with seismic 
sounding and electrical resistivity methods, the distance 
between the source and detecting elements is an independent 
variable. But, additionally, with the inductive coupling 
method, signal frequency is an independent variable. This 
additional degree of freedom makes it possible to collect 
a wealth of semi-redundant data that can be used to resolve 
experimental errors and uncertainties. And, because the 
SWEEP procedure promises to be fast and economical, one can 
make many traverses in different directions to distinguish 
between vertical and lateral variations in conductivity. 

2.1 Background 

Several seemingly unrelated developments and 
observations coalesced in a decision to undertake the 
feasibility study reported here. 

In the MIT Apollo 17 Surface Electrical Properties 
Experiment, conducted successfully on the moon in 
December 1972, electromagnetic energy was radiated from a 
small radio transmitter and horizontal electric dipole 
deployed at the landing site. A receiver mounted on the 
Lunar Roving Vehicle measured, and recorded on a portable 
magnetic tape recorder, the resulting field strengths as a 
function of distance. By interpreting the interference 
patterns produced by the various waves that travelled 
above and below the lunar surface, it was possible to 



determine some of the electrical characteristics of the 
lunar interior. One of the early objectives in the lunar 
experiment was to search for the presence of <n the 

subsurface of the moon. As a result of tho . measu? ments , 
it is now certain that at the Apollo 17 landing site, 
there is no free water at least to depths of several hundred 
meters . 

Prospecting for subsurface water also has increased 
urgency here on earth. We are aware of the serious problem 
caused by the continuing drought in the Sahel region of 
West Africa. This region, including the countries of Mali, 
Upper Volta, Senegal, Mauritania, and the southern reaches 
of the Sahara Desert, is now in the sixth year of subnormal 
rainfall that has caused great deprivation and widespread 
famine. In the hope of helping to alleviate some of the 
problems, we explored the possibility of adapting the lunar 
technology and experience to the purpose of finding subsurface 
water in arid regions of the earth. 

A direct transfer of the lunar technology is not 
feasible because, at the radio frequencies appropriate for 
an exploration of the lunar interior, the electrical losses 
in the earth are intolerably large. 

However, considerable attention has been given in recent 
years to the diffusion through the earth of electromagnetic 
energy at very low frequencies for Which the electrical 
losses can be very small. Several theoretical solutions for 
the fields produced by a loop in the presence of a layered 
conductive earth have been published. And experiments have 
been reported in which the conductivity of the subsurface 
was measured by observing the fields produced by a magnetic 
dipole at the surface. These experiments, for the most part, 
were handicapped by the difficulty of generating high power, 
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low frequency, signals with easily transportable 
apparatus . 

As part of another program, attention has been directed 
recently to the design of a variable voltage, variable 
frequency, power supply useful for converting energy stored 
in batteries to the form required to drive an electric 
motor in an electric vehicle. Coincidentally, a switched 
source synthesizer technique inspired by this purpose can be 
adapted readily to the frequencies and power levels 
appropriate for the inductive coupling method of measuring 
subsurface conductivity of the earth. 

So, the advent of a practical method for generating the 
required high power, variable frequency, source signals, 
the theoretical formulations describing the interaction of 
low frequency signals with the earth , the evident need for 
a prospecting technique permitting fast and economical 
survey of subsurface water resources in a large arid region 
of the earth, a considerable theoretical base result 
from a lunar experiment, and the practical implem.. ~on of 
the automated equipment used in the lunar experiment, 
combined to suggest that it would be worthwhile to examine 
the feasibility of an automated prospecting tool for use on 
earth . 

The results of this feasibility study, reported here, 
are most encouraging. 
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3.0 The Experiment Concept 

As with most remote sensing techniques, the fine detail 
of the rigorous analysis may be tedious and abstruse, but 
there is an intuitive approximation available that provides 
a powerful insight into the problem. In this section we 
describe the concept of the technique and give a physical 
explanation. 

For the mathematically inclined reader, this section 
should provide a useful basis for the more mathematical 
section on theory. For the reader who is not 
mathematically inclined, this section should provide a good 
basis for understanding the method and he may wish to omit 
the later theoretical section. 

Consider two coils in free space, such as shown 
schematically in Figure 3-1. One coil, the source, is 
energized with an alternating current. The other coil, 
the detector, is either coplanar with the source coil or it 
is perpendicular to the plane of the source coil and has its 
center o.. the plane defined by the source. With the source 
coil energized with an alternating current, there will be a 
time-varying magnetic field that can be detected by the 
other coil. With the configuration shown in Figure 3-1 (a), 
the detecting coil that is oriented in the plane containing 
the source coil ("horizontally") will be maximally coupled 
to the time-varying field while the coil that is oriented 
orthogonal to the plane ("vertically") will be completely 
decoupled and unaffected by the time-varying field. Thus, 
at the terminals of the horizontal coil there will be an 
alternating voltage having a magnitude dependent on the 
distance between the source and detecting coils and the 
magnetic moment, or intensity, of the source coil. But no 
such voltage will exist at the terminals of the vertical 


8 



Figure 3-1 (a) . A Source Coil With Two Detecting Coils in the Same 
Plane 


I 

I 

i 



Figure 3-1 (b) . A Source Coil With Two Detecting Coils Displaced 
From the Plane 
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coil because, in free space, this coil is decoupled from 
the time-varying magnetic field. 

If the detecting coils are displaced from the plane 
defined by the source coil, as in Figure 3-1 (b) , both coils 
will be coupled in varying degrees to the time-varying 
magnetic field and so an alternating voltage will appear at 
the terminals of both detecting coils. Calculation of the 
magnitude of these voltages is a straightforward geometrical 
problem for which the solution is well known. 

If the two coils are not in free space but are near a 
conducting medium, then the magnetic field will not be 
symmetrical about the source coil but will be spatially 
distorted by an amount that depends on the conductivity of 
the medium. The practical situation is that which occurs 
when the source coil is laid flat on the surface of the 
earth. See Figure 3-2. The field is distorted because 
currents are induced in the conducting material by the 
source coil and these currents modify the spatial 
distribution of the magnetic field. Due to this 
distortion, both detecting coils will yield an output 
voltage that depends on the nature of the medium on each 
side of the interfacial plane. 

The effect of this distortion can be visualized readily 
on the basis of a mathematically equivalent image of the 
source coil. The exact description of the image includes 
information on location, intensity, and time phase which 
are so adjusted that the superposition of the free-space 
field originating from the source coil and the hypothetical 
field arising from the image combine at the detector 
location to describe correctly the intensity, direction, 
and time dependence of the actual field. This image concept 
is illustrated in Figure 3-3 and the mathematical basis is 
described in Appendix D. 
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In actual practice, the situation is somewhat more 
complicated than the simple pictures of Figures 3-2 and 
3-3 might indicate because the presence of a conducting 
material not only affects the spatial distribution of the 
magnetic field but also affects the time phase of the field 
relative to that which would occur in a free-space 
environment. However, this can be accounted for by 
associating a time phase with the image coil so that the 
time dependent phase of the hypothetical field associated 
with the image source is different from the phase of the 
free-space field associated with the original source coil. 

The situation of interest to us is still more 
complicated because the earth is not a homogeneous body 
with one conductivity but generally consists of successive 
strata having different conductivities that are determined 
mainly by porosity and water salinity. Indeed, it is this 
layered structure and moisture dependent conductivity that 
we rely on to indicate the presence of subsurface water. 

Figure 3-4 illustrates the simplest, idealized, model. 
If the first subsurface layer has zero conductivity, that 
is, it is an infinitely resistive insulator, then no 
currents will exist in this layer. But if the second layer 
consists of water saturated rock that is highly conductive, 
that is, it is a relatively low resistance conductor, then 
the source magnetic field will induce substantial currents 
in the second layer. These currents, in turn, will create 
a magnetic field which, viewed from the surface, will 
appear to have originated from an image coil located at 
some complex depth. 

Although the image description of the magnetic field is 
a powerful conceptual tool, its application to the data 
analysis, as shown in Appendix D, is limited to large 
source-receiver separations. Since most of our data was 
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Figure 3-4 . The Image Concept for the Two Layered Case 





obtained at separations where the image formulation was 
invalid, we were forced to use a more exact integral 
formulation of the magnetic field. In the next section we 
snail describe the mathematical techniques used to formulate 
expressions for the magnetic field of a horizontal source 
loop on the surface of a multilayered earth. 



4.0 Theoretical Analysis 


In this section, we shall derive expressions for the 
magnetic fields induced by a horizontal current loop over a 

multilayer conducting earth in order to provide a 
theoretical basis for the analysis of the El Paso field 
measurements described in the subsequent chapters. 

We shall assume: 

a) a multilayered flat conductive earth, 

b) the permeability (y) of the earth is that of 
free space (yo), 

c) displacement currents can be ignored (i.e., 
the quasistatic approximation holds) , 

d) the horizontal current loop can be approximated 
as a vertical magnetic dipole of moment m 
(i.e., a point source), 

e) conductivity varies only with depth, 

and f) the current varies sinusoidally with time, with 

an implied e lwt time dependence in all expressions. 
The current source is placed on the surface of the earth (h=o) 
and at the origin of a cylindrical set of coordinates with 
our observation point P at (r,z) as illustrated in 
Figure 4-1. 


4 . 1 Magnetic Hertz Potential 

To simplify the analysis of electromagnetic boundary 
value problems, "gauge functions" are defined in terms of 
the measurable fields, the electric field vector E 
(volts/meter) and the magnetic field intensity vector H 
(amperes/meter) . The gauge we use for the current loop 
source is the magnetic Hertz vector n which is expressed 
in terms of E and H as: 


-% 
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H 


(V x V x II) 


(4 . 1A) 


E = iwy(V x n) (e 1U)t time dependence implied) (4. IB) 

The Hertz potential for a magnetic dipole is vertical. 
Therefore, the Helmholtz equation in layer 0 may be written 
in cylindrical coordinates as: 

V2n 0z + k2n 0z = “ m6(R) (4.2) 

where : 


6(R) = 2* j^oUrU dX (4.3) 

e 

In the lower (source free) layers, we can write* 

(V 2 + k. 2 )n. = 0 (4.4) 

1 1Z 

2 

where k. 2 = (io^uow - e^you )~ia^y 0 w (quasistatic 
approximation) and to determine the boundary conditions on 
H, equation (4.1) can be written in the form: 

- V 57W "z - H r> ■ 0 

-iwfl",- i E r = E 2 = 0 

The application of the normal boundary conditions for E 

and H on the Hertz vector n z yields the constraints that 

n and 311 /3z are continuous across the layers. We can 
z z 

now solve the Helmholtz equation in each layer and apply 
the boundary conditions to find II 2 at the observation 


(4.5A) 
(4 . 5B) 
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ooint P(r,z) above the surface of the earth. From the 
solution for II and equation (4.5) , we can readily deduce 
expressions for our measured field components H z and H r> 

The homogeneous solution for n z is of the form: 

( n z ) h = |j 0 (X r )« w - (Aa " fc * ) ** F (X) dX (4.6A) 

o 

where F(X) is an arbitrary function of X. 

The particular solution of the Helmholtz equation is 
of the form: 

ikR 

(II) = 22 (4.6B) 

2 p 4 ttR 


which can be put in the Sommer f eld integral form: 


'Vp * tBl 


ikR 


w 

V- l S e-°l z| J 0 (X r ) ax (4.7) 


where: 


2 2 4- 

a = (X - k ) 


Thus, the complete solution in layer 0, which contains the 
source and the observation point, is of the form: 


n 


0z 


TfT f J o (Xr) 



e -a°z 


+ f 0 ( A ) e 



dX 


(4.8) 


Each lower layer expression except the last must have 
downward and upward diffusing terms representing transmitted 
and reflected waves. 



4.2 Two Layer Case 


Consider now the two layer case where the lower layer 
extends to infinity. We can write expressions for II ^ and 
n 2z for layers 1 and 2 in the forms: 

00 

n iz = W~ f J 0 (Ar) [f!(^e aiz + f 2 U)e~ ai ^J dA (4.9) 

^22 = }Jo (Xr > dX u.io> 


Now we can apply the boundary conditions 

"oz - n iz ? 2 - 0 


n ss N @ z = ~d, = -d 
lz 2z ^ 1 


(4.11 Aj 
(4.11 B) 


and 


an 


Oz 


an 


lz 


3z 


az* 


@2 = 0 


(4.12 A) 


91 lz . 9 "2z 
az az 


@ z = -d 


(4.12 B) 


Note that the application of the boundary conditions affects 
only unknown terms f^(A). Thus, we can consider these terms 
independently of their integrals and can write: 
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at z = 0 


f n (X)- f. (X) - f,(X) = -/ 


ao 


f 0 (x > + st ; £ 1< X > - £ 2 (X) 


- V 


(4.13) 


ao 


and at z = -d 

f 1 (X)e“°‘ ld + f 2 (X)e +0tld = f 3 (X)e" a2d 


ai 


f 1 (X)e" aid - f 2 (X)e +aid 


(4.14) 


= a 2 f 3 (X)e 


-a 2 d 


To determine it Q z r the Hertz vector on and above the surface , 
we need to find fg(X) in terms of on, a 2/ and d. Solving 
(4.13) and (4.14) for fg(X) we determine that: 



(4.16) 
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Noting that 

2 2 r 4 - 

a 0 = (X ~ k Q ) = (X z - ia 0 yow t e 0 yow ) 

(X 2 + eoyotn 2 )^’ (4. 17 A) 

and ignoring displacement currents, we can write; 

CIO * x (4.17B) 

and (4.16) can be written in the j t.i: 



(4.18) 


Using equation (4.5) to find H z and H r we find that; 



(4.20) 


I 

I 

I 




Letting a = X in equation 4.7 we note that: 


ra * w~ J J o (Xr)e Xz dX 


(4.21) 


and 


3 

az 2 



= |X 2 J n (Xr)e’* Xz dX 


(4.22A) 


3 _ 

TrTz 


uo 

(^ = x2j l 


(Xr)e" Xz d 


(4.22B) 


Thus, we can write equation 4.19 in the form 


H = ~ m + 3m z + _5L 
Z 4ttR 3 4ttR 5 4ir 


00 

J? 0 (Xr)e- Xz A ! } 


Nai + 


(a i + X) (ai - a 2 )e“ 2aid - (on - X) (cu + aj) 


(ai + X) (a i + a 2 ) - (ai - X) (cti - a 2 )e 


and equation 4.20 in the form: 


^TTdj dX 

(4.23) 


„ _ +3mrz , m 

n r — - + 

4ttR 5 4tt 


[j 1 (Xr)e’ Xz X 2 


2 

✓ 


(a: + X)(ai - a 2 )e"’ 2ctld - (ai - X) (otj + a 2 ) 


(cxi + X) (a i + a 2 ) - (aj - X)(ai-'.a 2 )e 


r 5aTT 


d>. 


(4.24) 
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Following closely the analysis of Wait (1958) , we 
shall normalize equations (4.23) and (4.24) by dividing 
these expressions for H z and by the expression for the 
primary H z field (i.e., the first two terms of equation 
(4.23) and we shall, in addition, put these normalized 

expressions in forms convenient for numerical integration. 
Accordingly, the expression for the normalized magnetic 

field can be put in the form: 


h" , W*) 3 T.(A,B) + 1 (4.25) 

2 1 - 3(z/R>* 0 

and the expression for the normalized radial magnetic field 
can be put in the form: 


H n = 


(R/6) 


1 - 3 (z/R) 


• T. (A,B) - '-- rZ) 

2 1 - 3(z, 


(4.26) 


(z/R) 


where: 


T i (A,B) = Jr (D, g) g 2 e" gA J^gB) dg 


(4.27) 


where: 


is the Bessel function of the first kind of order i, 


and: 


R(D,g) = 1 - 


2g jug + e~ UD ) + V(1 - e" UD )l 


, / ti 2 . . _**\ / 1 ~ “UD x 


U (g + V) (1 + e ) + (U z + gV) (1 - e 
= 1 2q |(U + V) 4- (U - V)e “ UD ) 


(U + g) (U + V) - (U - g) (U - V) e 


■ = TTfr 


(4.28) 
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where: 


U = (g 2 + i2) 
V = (g 2 + i2 


. i 



r 2 = x 2 + y 2 
R 2 » r z + (z 


2 


02 

0i 



(4.29A) 



2 

) 

(4.29B) 

503 . 8 

(4.29C) 

(Cif) 2 

(4.29D) 


(4.29E) 


(4.29F) 


(4.29G) 

h) 2 

(4.29H) 


It should be noted that for all practical source 
receiver separations where z<<r, so that r-R* equations 
(4.25) and (4.26) reduce to: 


H z = (f)‘ V A - B) + 1 

and 

H ? -(f ) 1 t i < a ' b) - H 


(4.30) 


(4.31) 
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4.3 Integral calculation 


To evaluate equations (4.25) and (4.26) we have to 
evaluate the Bessel function integral of (4.27), For this 
purpose, we follow closely the numerical integration 
technique developed by Longman (1957) . 

First, we let: 


y = gB 


yA 

■ilA.B) = j <> B J i<y> d (£) 

o 

- -^7 R ( D '&) y 2 “P [~ y<2 r~^ ) J i 

O 

00 

= - F(y) J.(y) dy 
B 3 


(y) dy 


where : 


F (y) = r(d, J} y 2 exp (Z ^ - +h - } ~ 


I 


I 

I 

I 

I 


i 

I 

I 

I 

I 


The integrals T^(A,B) then are evaluated as the suit> of a 
set of integrals, each taken between successive zeros of 

J^y) * 


T. (A,B) = 


r 

z . , 
• i , 1 

?i,2 J 

< 

F (y) J i (y) dy + 

j 

F(y) J^y) dy + — > 

v < 

J 

z i,l J 


= - j F(y) J.(y) dy +2 




J 

0 


j=l 


F(y) J t (y) dy 


.1 r, 


Z if j 


J 


B ' 


S i,0<Y> + S S i,j (y) 
j=l 


(4.35) 


where: ' 

0, ( z • i ) , {?,. ~) are the successive values of y 

for which J\ (y) is zero 

and: z. . 

t 

(4.36A) 


S i>0 ^ ~ F(y) J^y) dy 


S i> j (y ) 


F (y) J i (y) dy 
z . . 


(4.36B) 
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To evaluate the integrals S . . (y) we use the 16 point 

1 i J 

Gaussian quadrature approximation: 


16 


(X) dX = w k f(X k } + m 


(4.37A) 


-1 


X=1 


or, transforming the limits: 
Q 16 


[ f (X) dX = W k ^ ^ Q - ~2 - P X k + Rm (4.37B) 


where and X^ are the Gaussian values listed in Table 4-1. 


So: 


s i f-'(y> 


Z i,jtl ~ Z i,j 


16 


2>k F( y i ,j,k ) J i (y i, j,k )+Rra 


k=l 


16 

S 

k=l 


= 2T c i,j.k F( yi,3,k ) + *” 


(4.38) 


where : 


C i,j,k - ~~ — W k 


v . ~ z i.j 

Jfi,j,k 2 


« , Z i-i+l ~ Z i,1 

X k + ^2 


( 4 . 3 9 A ) 


(4.39B) 
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Table 4-1 


Gaussian Integration Coefficients (W^) and Abcissas (X^) 


0.989400934991650 

0.944575023073233 

0.865631202387832 

0.755404408355003 

0.617876244402644 

0.458016777657227 

0.281603550779259 

0.095012509837637 

0.095012509837637 

0.281603550779259 

0.458016777657227 

0.617876244402644 

0.755404408355003 

0.865631202387832 

0.944575023073233 

0.989400934991650 


0.027152459411754 

0.062253523938648 

0.095158511682493 

0.124628971255534 

0.149595988816577 

0.169156519395002 

0.182603415044924 

0.189450610455068 

0.189450610455068 

0.182603415044924 

0.169156519395002 

0.149595988816577 

0.124628971255534 

0.095158511682493 

0.062253523938648 

0.027152459411754 
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We evaluated the expressions for the normalized magnetic 
field (equations 4.25 and 4.26) on an IBM S/370 Model 165 
computer using double precision arithmetic to evaluate the 
terms T^(A,B). This degree of precision (16 decimal digits) 
was found to be necessary to avoid roundoff errors because 
the evaluation hinged on small differences between large 
numbers; in retrospect, we realized that this problem could 
have been minimized by rephrasing the equations in terms of 
a half space solution and perturbations thereon. The 
integrals S. .(y) (equation 4.38) were summed over the first 

l 

100 zeros of the Bessel functions J^(y) and a one term Euler 
transformation was performed on the summation by subtracting 
half of the last term, ^Q Q (y). We determined early in 
the process that this brute force approach would yield more 
than adequate rigor for our purpose and, because of the 
pressure of time, we did not investigate the reduction in 
computer running time that we realized could have been 
achieved by reducing the number of zeros and using more 
terms in the Euler transformation. The term (1 - e UD ) 
was set equal to unity for values of UD greater than 39, 
because for this value Re|e” UD | 10**^ and, since the 

computer carries a double precision variable as a 16 
decimal digit number, values beyond this limit simply cause 
an underflow to occur without contributing anything useful 
to the results. 




4 . 4 Review 


In this chapter, we have derived expressions for the 
near surface normalized magnetic fields produced by a small 
current loop deployed on the surface of a layered conducting 
earth, and described the numerical integration techniques 
which we employed to evaluate the resulting expressions. 

In the next chapter, we describe field measurements 
made at two test sites near El Paso, Texas. Chapter 6 
then provides an analysis of these field measurements 
based on the preceding theoretical formulations. 

It should be noted that our data analysis to date has 
concentrated on a two layer theoretical formulation which 
proved adequate for the immediate purpose. As opportunity 
s rises in the future we would expect to incorporate into 
our computer program the multilayer analysis and inversion 
techniques explained in the appendices. 



5.0 A Field Trial 


A field trial was conducted near El Paso, Texas, during 
the month of October, 1973. This field trial was intended 
to be a rapid test of the concept and feasibility of the 
low frequency inductive technique for geophysical 
prospecting. For this reason, no attempt was made to 
refine the final hardware, but, rather, a collection of 
immediately available "black boxes" was assembled into a 
system with the sole criterion that the electrical 
performance had to be adequate to permit a realistic 
feasibility test. 

The El Paso area was selected for the initial field 
trial for several reasons. First, this area is quite 
arid and representative of areas where there is a real 
need for subsurface water prospecting techniques. Second, 
the area has been studied extensively in the past so that 
a comparative data base exists. Third, specific areas of 
interest are readily accessible. Fourth, and importantly, 
because the entire region relies on subsurface water, 
there are a very large number of operating wells that 
provide positive assurance of the depth and quality of the 
subsurface water. 

The climate is shown in Figure 5-1. The average annual 
rainfall in the area is only seven inches, a figure 
that is comparable with the rainfall in many areas of 
West Africa. The groundwater in the area has been described 
by Leggat, Lowry, and Hood (1963) . Although we used their 
study as general background information, the depths to water, 
as well as current conditions, were provided by Walter Meyer 
of the U.S. Geological Survey at the El Paso office. 
Additional geophysical data were available from Zhody 
(1969) . 
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(A) 



Average monthly temperature at 
El Paso, Texas 


(B) 
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monthly evaporation 
at Ysleta, Texas 


Precipitation range by month at El Paso, Texas 

Figure 5-1. Climate in vicinity of El Paso, Texas, 1878-1958. 

(From U.S. Weather Bureau records. 

After Leggat et al, 1963.) 
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In Figure 5-2 we show the location of the two sites 
that were used: in the Messila Valley at the City 

Well Field and on the Santa Teresa Mesa. These locations 

were chosen because they offered good calibration of 
depths to water on the basis of numerous wells, and 
because they offered markedly different depths to the 
water table. The site used in the Messila Valley had a 
depth to the subsurface water of 8.6 meters (28 feet), 
and the Santa Teresa Mesa site had a subsurface water depth 
of 100 meters (327 feet). 

For the October field trial, the 'source’ coil 
consisted of 12 turns of #10 wire laid on the ground in a 
square with a side of 25 meters, and energized with 
currents ranging from 15 amperes at 60 hertz to 
5 milliamperes at 20 kilohertz. The ’detector' coil 
consisted of 1200 turns of #26 wire wound on a one meter 
square wooden supporting structure that could be oriented 
in either a horizontal or vertical position. The 
magnitude and phase of the observed signal was measured 
by means of a commercial phase-locked- loop synchronous 
detector using a reference signal coupled from the source 
to the detector by a shielded, twisted pair, cable. 

At both sites , data were collected with source-detector 
separations ranging from 20 to 500 meters. A photograph 
of the apparatus is shown as Figure 5-3. 

Figure 5-4 shews the soiree energizing equipment which 
included three motor driven alternators and an 
oscillator-amplifier combination mounted in a pickup truck. 
One 3.5 kVA gasoline driven alternator served as a portable 
power supply for other equipment and also provided a 
current of 15 amperes at a frequency of 60 hertz that was 
coupled directly to the source coil. A second, 1.2 kVA, 
gasolina driven alternator provided a current of 
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3.1 amperes at a frequency of 400 hertz that was coupled 
directly to the source coil, and a third gasoline driven 
alternator provided a currr- 4 - of 1.2 amperes at a 
frequency of 800 hertz that was similarly coupled directly 
to the source coil. Other frequencies of 1, 2.5, 7, 16, 
and 20 kilohertz were derived by using an oscillator to 
provide input to a high power commercial audio amplifier; 
currents supplied to the source coil ranged from 800 
milliamperes at 1 kilohertz to 5 milliamperes at 
20 kilohertz. 

The detecting apparatus is shown as Figure 5-5. A 
wooden framework supported a coil consisting of 1200 turns 
of #26 wire wound in four separated layers of 300 turns 
each on a square wooden form with a maximum cross-sectional 
area of one square meter. The coil structure was mounted 
on pivots so that it could be oriented in either a 
horizontal or a vertical plane; final adjustment was 
eccomplished with the aid of levelling bubbles. At the 
higher frequencies, data were measured using only the 
inside layer of 300 turns; this layer had a self resonant 
frequency of 35 kilohertz and a Q of 37. At ranges 
greater than 300 meters, where the higher frequencies wore 
no longer of interest (because the signal level was too 
low for measurements) , measurements were made only at 
frequencies of 1000 hertz or less: in this regime we 

used the full 1200 turns which had a self resonant 
frequency of 6 kilohertz and a Q of 17 . The battery 
operated synchronous detector was mounted on the framework, 
in a position chosen to be symmetrical with respect to the 
coil, and a phase reference signal was provided by a 
shielded, twisted pair, cable extended between tht_ source 
and the detector apparatus. 
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Figure 5-3. Photograph of apparatus used in field trial. 
The source energi 2 ±ng equipment was 
mounted in the tarpaulin shaded, pickup truck. 
Part of the horizontal source coil can be 
seen in the foreground, and the detector 
apparatus can be seen in the midgrourd beyond 
the hood of the pickup truck. 




Figure 5-4. Photograph of source energizing equipment. 

Three alternators (60, 400, and 800 Hertz) 
were carried in the truck w*th an oscillator , 
power amplifier, oscilloscope and switching 
box mounted on the tailgate. 




Photograph of detector apparatus. 

1200 turns of #26 wire were wound on a 
wooden form that could be rotated to either a 
vertical or horizontal position. The phase 
sensitive detector was mounted symmetrically 
with respect to the coil. A phase 
reference was provided by a shielded, twisted 
pair, cable extending from the source 
apparatus. Small hydraulic jacks were used 
to level the apparatus at each recording 
station. 
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Measurements were made at logarithmically spaced 

stations (10 1 * 3 , 10 1 * 4 , 10 1 * 5 , .*.) from 20 meters ( 10 1 * 3 ) 

2 7 

to 501 meters (10 * ) as the detector apparatus was carried 
in a straight line away from the source. At each station 
the magnitude and phase of the voltage existing at the 
terminals of the coil were recorded, first with the coil 
vertical and then with the coil horizontal. 

Six frequencies were used at each of the shorter range 
stations but as the range increased the signal-to-noise ratios 
of the higher frequencies progressively decreased below a 
useable level. Accordingly, the number of frequencies was 
reduced until, at the longer ranges, data were recorded only 
at 60 and 400 hertz. However, since most of the information 
at long ranges is related to deep layers and can be gleaned 
from the low frequency data, the absence of higher frequency 
data at the longer ranges was not detrimental to the 
feasibility test. 

It should be noted that the purpose of the field trial 
was to test the concept and feasibility of the technique 
for geophysical probing. The apparatus only vaguely (if at 
all) resembled a final operational system. We needed about 
three hours to deploy and recover the source coil; while 
the detector apparatus, which weighed more than 100 pounds, 
had to be manhandled with considerable difficulty. During 
this field trial, one data run consisting of fifteen 
stations ranging from 20 to 500 meters from the source 
occupied two full days of effort with an average time of 
16 minutes for each elementary datum point (one frequency 
at one station). In an operational system, of course, 
smaller and more efficient source and detector coils would 
be mounted permanently on vehicles capable of traversing 
the type of terrain of likely interest. With proper 
instrumentation of this type we expect to obtain a complete 
data run in one hour. 
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6.0 Data Analysis 

For our feasibility trial field measurements near 
El Paso, Texas, we chose sites of known geology. 
Specifically, we chose sites where the depths to the 
water table were defined clearly by existing wells, and 
where we could work over large, apparently homogeneous, 
areas which appeared to be uniform in all directions from 
the source coil. Thus, in the data analysis we could 
fairly assume negligible lateral variations in 
conductivity, i.e., a (r ,z) =c(z) . 

Although we realize that variations in conductivity 
with depth are rarely abrupt in the type of soils of 
interest in this investigation, we chose to assume in the 
data analysis that the subsurface of the test sites 
consisted of discrete, homogeneous, layers of varying 
conductivities. That is, we chose to assume in the data 
analysis that the conductivity with depth could be 
approximated adequately as a discretely varying parameter, 
a(z) = a^(z) where i is the layer number. This 
assumption of a layered medium greatly simplifies the 
theoretical problem, allowing us to use the analysis of 
Chapter 4 and Appendix A for the calculation of the 
fields induced by an oscillating magnetic dipole near a 
multi-layered conducting earth.. 

We assumed in our theoretical models that the 
permeability of the subsurface was the same as that of 
free space, y o , since the permeability of moist rocks and 
soils is very close to that of free space except when large 
quantities of ferrous minerals are present. And, finally, 
we assumed that displacement currents could be ignored 
for our low source frequencies. A summary of the above 
assumptions is listed in Chapter 4 . 
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6.1 Data Quality and Format 

During our field tests, we obtained measures of the 

absolute magnitude and phase (with respect to the phase of 

the source) of the radial and vertical components of the 

magnetic field, H and H , for the frequencies and ranges 

i z 

described in Chapter 5. At each frequency we took data at 
increasing separations until limited by noise in the 
receiver. Data points at the maximum range for each 
frequency, therefore, can be considered to be of suspect 
quality because they were derived from measurements made 
with a very low signal-to-noise ratio. 

In our analysis we matched theoretical curves 
derived from expressions for the normalized field 
components (see Chapter 4) with plots of the measured 
values. The quantity actually measured in the field is 
the voltage induced in the receiving loop when tne source 
loop is energized. 







For ranges much greater than the height of the receiving 
coil (the usual situation of interest) the measured voltage 
is related to the normalized components of the magnetic 
field by the equation: 



where i 


iam g IN lAiN^Ag 
4»r 3 



(- 1 ) 


1/2 


( 6 . 1 ) 


w = 



source frequency in radians/sec. 

-7 

permeability of free space = 4txl0 

magnitude of source current in amperes 

number of turns of source loop 

2 

area of source loop xn meters 

number of turns of receiving loop 

2 

area of receiving loop in meters 

normalized component of the magnetic field in 
amperes/meter 


r = separation of source and receiver in meters. 
Normalisation of the data, therefore, is achieved by 
multiplying the measured voltages by the factor: 

4-irr 3 

iu)HoINiAiN2A2 


We should note here that, for convenience in the field 
trials, we arbitrarily chose the orientation of the 
receiving loops and the terminal connections so that all 
phase measurements approached zero degree at short ranges. 
Accordingly, for consistency in comparing the data with the 
corresponding theoretical expressions we shifted the phase 
of the theoretical expression for the radial magnetic field 
component by 180 degrees; this is simply equivalent to 
rotating the vertical receiving loop by 180° or reversing 
the terminal connections to the loop. 
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6.2 Analytical Development 

Initially, we based our analysis on the approximation 
that the subsurface of the test sites could be considered 
as consisting of only two layers of different 
conductivities, with the second layer extending to depths 
below the detection capability of our apparatus. As we 
shall show, this approximation proved adequate for the City 
Well Field site which had a water table at a depth of 
9 meters, but inadequate for the Santa Teresa Mesa site 
which had a water table at 1QQ meters; for this latter 
site f we were forced to a three layer model. Time and 
budget considerations prevented us from extending the 
theory and the computer programs to the multilayer case 
in the manner described in the appendices. Instead, we 
chose to handle the analysis of the Santa Teresa Mesa 
data by using two different two- layer models in different 
frequency regimes, one model defined by the first and 
second layers, and one defined by the second and third 
layers. It might be noted here that the possibility of 
separating the problem into two (or more) parts in this 
manner is one of the attractive features of the 
inductive coupling technique. 

We started the analysis using the simple image theory 
outlined in Appendix D as the basis of our calculations of 
the normalized field components H and H but it soon became 

XT Z 

evident that the image concept could be applied accurately 
only at ranges large with respect to the effective skin depth 
of the upper layers (see Appendix D) and that this restriction 
limited the usefulness of this approach. 
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We turned then to the integral solution described in 
Chapter 4, in order to provide more accurate theoretical 
curves. The integral solution, indeed, provided a far 
more accurate description of the data than did the image 
solution, especially for separations less than the effective 
skin depth of the upper layers. 

However, after laborious attempts to match theoretical 
curves with data, we found that we still could not achieve 
plausible fits, and at that point we reexamined the subject 
of the field measurements and normalization. We determined 
eventually that our initial analysis problems were due 
largely to: 

a) bias in the measured phase due to a phase error 
in the reference signal; 

b) phase and amplitude errors in the field 
measurements caused by operating at frequencies 
too close to the resonant frequency of the 
receiving coil; 

c) errors in the normalization process of converting 
raw data into a convenient form for analysis, due 
primarily to an imprecise measurement of the source 
energizing current and the areas of the loops; 

and d) near field errors due to the finite size of the 
source loop . 


44 


All of these problems, save possibly d) , were circumvented 

by simply taking ratios of the magnitudes of the vertical 

and horizontal field components and differences between 

their individual nhase measurements. By considering 

ratios of the magnitudes of the field components, most of 

the errors in the normalization process disappear and only 

errors in the actual knowledge of range, and short term 

variations in the amplitude of the source energizing 

current, remain to hamper matching the theoretical and 

data curves. Additionally, the effects of receiver loop 

tuning, if only a function of frequency (not position) , 

disappear. Finally, by considering phase differences 

(4> -<}>), errors due to phase shifts in the 

reference signal are eliminated, since we can write the 

measured phases in terms of the real phase plus a bias 

error which is the same in both the radial (<f> r ) and 

vertical (<(> ) phase measurements, 
z 

Thus, our final choice of data format is bilogarithmic 
plots of the magnitude ratio | H /H | versus separation for 
a discrete set of frequencies, and semilogarithmic plots 
of the phase difference (<j> - 4> ) versus separation for a 

Z L 

discrete set of frequencies. Only near field errors have 
been neglected in the analysis but these errors occur 
only at combinations of short ranges and low frequencies 
where little useful information can be inferred from the 
data. 

Before presenting the results of our analysis, we shall 
first present the technique we ultimately used to 
determine the nature and bounds of the variation of 
conductivity with depth so that we could judiciously 
choose an initial conductivity model for the El Paso 
field sites. 
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6,3 An initial Interpretation Scheme 

One method of presenting the data in a form diagnostic 
of conductivity variations with depth is to plot, on a 
bilogarithmic scale, the magnitude of the ratio H /H 

Z IT 

versus the quantity; 

(r2 ™‘ £,V2 "fej * e£fl/2 

For, if the earth could be considered as a homogeneous 
half space characterized by a single conductivity a, then 
a e ££=CT and plots of data, for various source-receiver 
separations (r) and frequencies (f^) , would fall on a 
single curve. Conversely, if the conductivity were a 
function of depth, plots of the data at each range for a 
fixed set of frequencies would tend to fall on separate 
curves, the shapes and positions of which would be 
dependent on the nature of the conductivity variations. 

For example, if the conductivity increased monotonically 
as a function of depth, cr ^ would increase with 
increasing r for a fixed set of frequencies and, 
correspondingly, the curves would be translated to the left 
for increasing values of r. Conversely, monotonically 
decreasing values of conductivity with depth would result in 
the ratio curves being translated to the right with 
increasing values of separation r. 

For the two layer case where the upper layer is less 
conductive than the lower layer, the ratio data curves, 
therefore, will tend to merge, as the separation is 
increased, on a single curve characteristic of a half 
space whose conductivity is that of the lower layer; 
whereas, as the separation is decreased, the ratio data 
curves will tend to merge with a half space curve 
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characterized by the conductivity of the upper layer. 

For the multilayered case, multiple half space 
asymptotes will occur from which an estimate of the 
vertical conductivity profile can be determined. 

We realize that an analogous technique could be used to 
derive an initial interpretation based on the phase data; 
however, time and budget considerations prevented us from 
pursuing this possibility. 

We conducted our analysis, therefore, by making an 
initial assessment of the vertical conductivity pi ■'file 
based on magnitude ratio alone , and we then obtained a 
final corroboration of the interpretation by comparing the 
data curves for both magnitude and phase with model curves 
generated on the basis of these initial estimates. 



6.4 City Well Field Results 


Applying the inspection technique outlined in the 
preceding section, we can infer that the average value of 
conductivity tends to increase with depth (Figures 6-1 and 
6-2) . We note in the Well Field data that the ratio curves 
tend to merge on a single curve for ranges greater than 
150 meters to the maximum range of 501 meters. By 
matching a half space model curve to the data for ranges 
greater than 150 meters, we can readily infer a 
conductivity^ of 0.1 siemens/meter. Thus, we can assume 
that a fairly thick layer of 0.1 siemens/meter must exist 
at a depth less than 150 meters. As range is decreased 
below 150 meters, each subsequent data curve tends to move 
closer toward a single curve representing a conductivity of 
the order of 0.01 siemens/meter . Additionally, we can deduce 
that the thickness of the upper layer must be less than 25 
meters since the data curves have not merged on a single 
curve for ranges as small as 25 meters and the presence of 
the lower layer is still noticed in the data. 

To arrive at our final model of conductivity as a 
function of depth, we then plotted two layer theoretical 
model curves in the forms of magnitude ratio l H z / H r l versus 
separation on a bilogarithmic scale and phase difference 
(<ji - $ ) versus separation (r) on a semi-logarithmic 

Z IT 

scale. We chose model parameters consistent with the 
asymptotic behavior of the data curves, as described in 
the previous paragraph, and compared our theoretical 
models with the data plotted i:i the same forms. 


TT) The term "siemens” and the symbol "S" have been adopted 
for the unit of conductance previously referred to as the 
"mho". Use of this term has been approved by the IEEE, IEC, 
ISO, and Conference Generale des Poids et Mesures. 
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Normalized Range ( g - g ^ y ) 
Figure 6-1. City Well Field Data - I 
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After performing several parameter iterations, we 
deduced the two layer model illustrated in Figure 6-3a 
though wo realized that the upper layer conductivity of 
0.01 siemens/meter is only an average value for the upper 
9 meters and that the conductivity cf the upper few meters 
is a much lower value. 

The curve matches on which this mocel is based are 
present jd as Figures 6-4 through 6-9, and as we can see in 
the high frequency curve matches (6-8 and 6-9) , the best 
fits occur for a marked gradation in conductivity for the upper 
nine meters. (Note that 6-8 and 6-9 use depths of 6 and 4 
meters for the upper layer.) Our interpretation of this 
vertical conductivity profile is that the rapid increase in 
conductivity in the upper nine meters is due to an 
increase in moisture with depth terminating in the water 
table at nine meters where no further increases in 
conductivity are detected. We also infer from the high 
level of conductivity below nine meters (0.1 siemens/meter) , 
that the water tends to be brackish. 

For comparison with our interpretation of the City Well 

Field site, we present as Figure 6-3b a smoothed set of 
well log resistivity measurements taken near our test site 

and supplied by Walt Meyer of the U.S.G.S. in El Paso, Texas. 

Clearly our model is consistent with the depth to the 
water table and the average value of conductivity of the 
layer of soil directly above the water table. Our 
assessment of the conductivity of the lower water 
saturated lo.yer appears consistent with the well log to 
approximately 70 meters beyond which the U.S.G.S. data 
indicates that the conductivity tends to decrease with 
further increases in depth. Thus, we must conclude that we 
could not detect the decrease in conductivity with the 
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frequencies and separations we used for depth sounding 
the City Well Field site. Wait and Fuller (1972) 
have described analytically this limitation inherent in the 
induction sounding technique when attempting to detect 
deep resistive layers overlain by more conductive layers. 
However, although we were unable with our feasibility test 
instrumentation to detect the presence of the resistive 
zone below 70 meters, we were able to detect the 
conductivity contrast attributable to the presence of 
ground water. 
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Figure 6-8. Model Curves Compared with City Well Field Data: 
7000 Hertz. 
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Figure 6-9. Model Curves Compared with City We]i Field Data: 
20000 Hertz. 
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6.5 Santa, Teresa Mesa Results 

Applying the inspection technique outlined in 
Section 6.2 to the Santa Teresa Mesa data, we can readily 
infer that the average value of conductivity at the mesa 
site tends to increase with depth (Figures 6-10 and 6-11) » 

We note in the Mesa data that the ratio curves between the 
separations of 100 to 200 meters tend to merge on a single 
curve from which a conductivity of 0.05 can be inferred by 
matching a half space model to the data curves. Thus, we 
can assume that a layer of conductivity 0.05 siemens/meter 
exists at a depth much less than 100 meters. For ranges 
less than 100 meters, the data curves, though ill defined, 
seem to approach a conductivity of 0.01 siemens/meter or less 
and we can deduce that the layer thickness must be less 
than 25 meters, since the upper layer half space asymptote 
is not even approached at the separation of 25 meters and 
the frequencies of 7000 and 16000 hertz. Additionally, we 
can see that for separations greater than 200 meters, the 
ratio data curves again tend to translate to the left 
indicating a third layer of higher conductivity below the 
0.05 siemens/meter layer. 

Matching half space curves to the data for large 
separations, we infer a third layer conductivity of 
0.1 siemens/meter or larger. We also infer that the third 
layer is at a depth of the order of 100 meters. We can 
make this inference because, for a conductivity of 0.05, 
the skin depth in the second layer is about 70 meters at 
400 hertz and about 170 meters at 60 hertz (the two lowest 
tes^: frequencies) and it would not be possible to sense 
the third layer at all at depths materially greater than 
the skin depth in the second layer. 
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We again plotted our theoretical models in the forms of 
magnitude ratio and phase difference versus range using the 
parameter values deduced with our initial inspection 
technique and compared our model curves with the data. 

We found that only 60 hertz was diagnostic of the 
depth to the third layer whereas the higher frequencies 
were most diagnostic of the upper layer. After performing 
several parameter iterations , we deduced the three layer 
model illustrated in Figure 6-12. The curve matches on 
which this model is based are presented as Figures 6-13 
through 6-18. We interpret the upper thin layer to 
represent fairly dry soil, the middle layer to represent 
moist but not saturated soil , and the lower layer to 
represent saturated soil or the top of the water table. 

Our interpretation compares favorably with data from nearby 
wells for which depths of 93 to 102 meters to the water 
table are reported by Walt Meyer of the U.S.G.S. in 
El Paso. Additionally, the conductivity values we 
deduced for this site are consistent with the average 
conductivity values for the upper 100 meters as reported 
by Zhody (1969) for test sites in the El Paso area. 
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Figure 6-16. Model Curves Compared with Santa Teresa Mesa Data: 
2500 Hertz. 
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Figure 6-18. Model Curves Compared with Santa Teresa Mesa Data 
16000 Hertz. 
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6.6 Curve Matching Sensitivity 

To clearly illustrate how well our model curves fit 
the data curves, we now present, as Figure 6-19 through 
6-21, overlain model curves representing the effects of 
perturbing the parameters o 1# o 2 , and d individually at 
a source frequency of 1000 Hertz. We should note that 
we used this iterative process of varying one parameter 
at a time throughout our study in order to obtain the 
best possible fit to the data within the limits of our 
initial assumptions . 
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Figure 6-20. Model Magnitude Ratio and Phase Difference versus 
Range at 1000 Hertz, with Variations in the 
Conductivity of the Upper Layer ( 01 ) . 
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7.0 Recommendations 


There seems now to be little doubt that the inductive 
coupling technique can be developed into a practical system 
for the quick and definitive determination of subsurface 
conductivities at depths extending to a few hundred meters. 
And, because subsurface water in arid regions usually 
results in sharp conductivity contrasts that are 
particularly easy to detect, this technique shows promise 
as a useful tool for mapping large regions in order to 
identify specific places from which subsurface water might 
be recovered economically. 

Work which needs to be done to bring this system to 
fruition includes extension of the theory to 
nonhomogeneous situations, the development of rapid 
interpretation procedures suitable for real time use in the 
field, the fabrication of automated operational hardware, 
and field tests designed to test the system thoroughly 
in a realistic environment. These various items are 
discussed below. 

7 . 1 Theoretical Analysis 

We chose to conduct the initial feasibility test in a 
geologically simple area where the subsurface structure 
was known to be homogeneous and isotropic, and to have 
only two or three layers. This made possible several 
simplifying assumptions in the theoretical work which, 
consequently, may not be applicable to all areas of 
potential interest. 

The theoretical work should be continued, therefore, 
in order to generalize the theory to more complicated 
geological structures having many layers that are not 
uniform, and to structures better characterized by 
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continuous qradations of conductivity rather than 
discrete layers. 

Additionally, in our data interpretation and analysis 
we used a large computer facility to make lengthy (and 
expensive) numerical integrations in order to generate 
specific model curves which we then compared with plots of 
the experimental data- This forward, curve matching, 
procedure initially had great advantage because it 
illuminated some discrepancies and biasses in the data, it 
directed our efforts toward the selection of reasonable 
models, and it enabled us to make judgements about the 
validity of some approximations and about the rigor with 
which numerical computations had to be performed. 

Howe'er, using a large computer to generate specific 
model curves which are later compared with the experimental 
data obviously is not a long term solution to the problem. 
The procedure is cumbersome and does not make efficient 
use of the computer capabilities; it requires a great deal 
of intuitive, and subjective, manual intervention; and it is 
clearly not suitable for real time data interpretation in 
the field. 

Work needs to be done in this area, therefore, to 
develop direct data inversion techniques (as outlined in 
Appendix C) , to simplify the numerical calculations (as 
outlined in Appendix B) , to generate catalogs of curves 
and to ide tify specific features of the data yielding 
rapid approximations to the subsurface conditions. And, 
importantly, this additional theoretical work is needed to 
help define system operating procedures that maximize the 
amount of interpretable data collected with a given 
expenditure of time and effort in the field. 


7.2 Operational Procedure 

The ultimate SWEEP system will employ two self-propelled 
vehicles, one carrying the source and one carrying the 
detector, which can move independently through a region 
that is to be mapped. 

With both vehicles stopped, data van be collected in a 
form convenient for immediate interpretation by sweeping 
the signal frequency or by\arying the frequency in discrete 
steps. Additionally, data can be collected as the 
separation between the vehicles is changed or as both 
vehicles move simultaneously along, or parallel to, a 
traverse line. 

With this many degrees of freedom, there obviously is 
opportunity for considerable ingenuity in devising an 
optimum operational routine that permits the rapid 
identification of areas of likely interest and the 
subsequent checking of postulated models. 

7 . 3 Equipment 

In previous attempts to exploit the inductiv*, coupling 
prospecting technique, one of the most difficult hardware 
problems, and, indeed, the one which usually has been an 
intolerable handicap, has been the implementation of a 
high power, variable frequency, source power supply. To 
provide a large enough system signal-to-noise ratio to 
permit measurements to depths of the order of a few 
hundred meters it is necessary to drive the source coil 
with currents in the order of 50 amperes at power levels 
in the order of 15 kilowatts at frequencies ranging down 
to about 10 hertz. 

In Appendix E, we describe a technique for solving 
this problem. A switched source synthesizer, developed 


for another purpose, makes possible a very efficient 
DC-to-AC inverter that creates a quasisinusoid by 
connecting storage batteries in series-parallel 
combinations. Since the full system power for SWEEP is 
required only in short spurts, just long enough to obtain a 
datum point, the storage batteries can be recharged with no 
difficulty from a modest sized, continuously running, prime 
power source such as a gasoline d 'iven alternator. The 
switched source synthesizer technique described in 
Appendix E, therefore, makes the small, high power, low 
frequency, portable power source a practical reality. 

As part of an operational system, we shall need a 
coherent detector of a form similar to that used in the 
initial feasibility trials, and a telemetry link to ensure 
that the source and detecting vehicles operate in unison. 
These items, however, require no technological advances 
but can be implemented readily with existing hardware 
techniques . 

7.4 Field Trials 

To date, we have obtained field data only from two 
sites near El Paso, Texas. Tnis experience should be 
extended as soon as possible, preferably by progressing 
from areas with well known geology to surveys of the 
unknown. 

Experience during the initial feasibility test showed 
(once again) that a great deal of time can be wasted 
fixing engineering problems in the field. For this reason 
the field trials should be conducted with a well engineered, 
reliable, system that has been thoroughly tested pi ’or to 
the actual field trials. Presently we think that we 
should instrument a pair of vehicles, one for the source 
and one for the receiver, and perform three local 
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shakedown engineering tests before embarking on full 
scale field trials elsewhere. 
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Appendix A; Multilayer Analysis 

In chapter 4 we described the analysis and calculation 
of the magnetic field induced by a magnetic dipole over a 
two layer earth. The mathematical procedures of Chapter 4 
of matching boundary conditions for the Hertz vector and 
solving simultaneous equations for the reflection and 
transmission coefficients could be extended to three or 
more layers. However, a simpler approach is to put the 
wave equations into the form of transmission line equations, 
considering each layer as a cascaded segment of a 
transmission line. 

Ward (1967) applied the transmission line analogy to a 
loop over a multilayer medium and his results are listed 
below. 
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As in equations (4.19) and (4.20) we can write 
expressions fo- H 2 and R for a current loop on a 
multilayered conducting earth, as shown in Figure 4-1, 
in the forms: 


H = 2- 
z 4tt 


R(X) X 2 J n (Xr) e" Xz dX 


IA.1) 


and: 


H r = I? j R(X> X 2 J i (Xr) e " XZ dX (e +lwt implied) (A. 2) 


\ here: 


R(X) + 1 = k ° - y - (A. 3) 

K 0 *0 

Z + K. tanh a.d. 

Z = K . -i i iJL (A. 4a) 

K 1 + Z 1 tanh a l d l 


Z 

n-1 



Z + K tanh ad 
n n n n 

K + 2 tanh ad 
n n n n 


(A. 4b) 


Z 

i-1 



(A. 5) 


where is the characteristic impedance of each layer 
and is represented by the equation; 



(A. 6) 


where : 


I 


Thus we have presented expressions for the magnetic 
field components induced by a magnetic dipole (current 
loop) over a multilayered earth. These expressions can 
be solved numerically using the techniques of section 4.3. 
However, for future work we propose the use of the simpler, 
less costly, numerical technique described in Appendix B. 



Appendix B; Integral Approximation 


We recall from Appendix A that the multilayer 
integral forms for the magnetic field H z and H r can be 
written as: 

00 

H z = 3F" J X * J 0 (Xr) e ~ XZ R(X) dX (B,1) 

O 

and: 

00 

H r = W“ | X2 J l (Xr) e " XZ R(X) dX (B.2) 

O 

Madden (1974) has simplified the numerical integration of 
equations of these forms by defining the kernel R(A) in the 
series form: 


R(X) 



(B.3) 


where b^ are a preselected set of coefficients the 
specification of which is based on the expected range of 
values of X for a given source receiver separation and 
conductivity-layering model. 

Now, we can write H z (or similarly H r ) in the form: 

CO 

3 2 -b . X , 

H = — - l e 1 a. J n (Xr) e Az dX (B.4) 

Z 3z 2 i 10 

« 

O 
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Reversing tne order 


H 


z 


00 

r 



o 


of integration and summation we have: 
~ (b . -tz ) X 

e 1 J Q (Xr) dX 


A 


- I a. L 
i 


3 b i 2 


[ (r 2 +b i 2 ) z 


<r 2+b i 2 )tJ 


(B.5) 


and if we define the elements of the vector a as and 
the vector c as c i where: 




tr 2 +b. 


2 )7 


(B.6) 


we can write: 


H 


z 


T 


a c 


(B.7) 


The next step, then, is to determine the vector a. 
Recalling, from B.3, that: 

-b.X 

R(X) = Z e 1 a. 

i 1 

and making the approximation that the continuous function 
R(X> can be represented by a discrete set of R(X^) which 
span the necessary range of values of X for an accurate 
approximation to the integral, we can write: 


and in matrix notations 


Rj = 


or 


R = Ba 


(B. 9) 


To find the vector a we premultiply both sides of 

T 

equation (B.9) by B=(B*) such that; 


&R = &Ba (B.10) 

and: 

a = (&B)' 1 &R (B.ll) 

where : 

B* is the conjugate of B 
B = (B*) is the transpose of B* 

and the use of & insures a unique matrix inversion to 
find a. 

Accordingly we can write; 

H z = [(&B)" 1 &r] T c (B ,12) 

Thus, we have a least squares approximation to the solution 
of the forward problem which can be used to reduce 
computation time drastically. An additional attraction to 
this shortcut integration approach is that much of the 
computation can be reused for the calculations of the data 
derivatives necessary for the implementation of the 
inversion scheme described in Appendix C. 


Appendix C; Data Inversion 

Data inversion is an iterative process in which 
perturbations in the_ parameters of a theoretical model are 
used to minimize, in the least squares sense, differences 
between the calculated model values and thw actual data. 

If the theoretical model closely approximates the 
physical situation, we can write: 


9tt z 

H (measured) = H (model) }+ (p. - p (model)) (C.l) 

2 3Pi 1 

where (p i - p(model)) are the variations from ths original 
model necessary to fulfill the equality of (C.l), or in 
matrix form: 


* H z 

AH = H (measured) - H (model) - * — Ap 
z z dp^ r 


(C . 2) 


and in the equal parameter weighting form: 


Jln-i 


H (measured) 
z 


H (model) 

t z 


a. 


3 (Jin H (model; ) 


J 3 Un p^ 


Lin p^ 


(C . 3) 


where p^ represents the i unknown parameter values of the 
model. Consequently, we have a set of equations in the 
matrix form 


Ay = d 


where y represents the variations in the parameters 

d the difference between observed and calculated H 

z 

and A the general m by n sensitivity matrix of the 
partial derivatives of the model equations with 
respect to the unknown parameters. 
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We wish to find y for which d goes to zero or a 
minimum. Thus, wc need to find the generalized inverse of 
A in terms of a set of nrn-zero eigenvalues. 

Following Professor Madden* s course notes (1974) , we 
find the generalized inverse to be of the form: 

y = (AA + £*E)" 1 A d (C.4) 

where E is the unit matrix and c is an arbitrary bias term 
necessary to prevent overemphasis of near zero 
eigenvalues of &A. 

Briefly, the above inverse formulation circumvents 
inversion difficulties due to small or zero eigenvalues 
ana squelches non-unique and incompatible solutions. The 
determination and iteration of the eigenvalues of the 
matrix A is the key factor in the realization of data 
inversion, since as we change y to minimize d, the matrix 
A also changes. We shall show now that much of the 
calculation of the sensitivity matrix has already been 
accomplished in the solution to the forward problem as 
calculated in Appendix B. 

C.l Sensitivity Matrix Calculations 

Recalling from Appendix A that the vertical magnetic 
represented as: 

J Q (Xr) X 2 e" Xz dX (C.5) 


field H z can be 




R(X) 


90 




m 


m 






and noting that all -he information for the partials of 
H z with respect to the unknown parameters and d^ is in 
the kernel R ( X ) , ve can write: 


J n ( Ar) A 2 e" Xz dA (C.6) 

3 Pi 

where p^ represents the unknown parameters a. and d^. 
Recalling from Appendix B that: 

H z = a T c (C.7) 

we can write the approximate sensitivity matrix as: 

— - - -2_ <a T c) - /(&B)' 1 S II l c (C.8) 

# PI spi t t) 

Thus, for each forward problem solution we can readily 
determine the corresponding sensitivity matrix, if we can 
analytically determine the partials of R(A) with respect 
to the parameters d^ and 

Instead of finding the partials of R(A) directly, we 
shall find the partials of 1+R(>) which is a simpler form 
to differentiate. We shall show that the partial 
differentials can be put in a simple recursive form so that 
the partial derivative with respect to any layer thickness 
or conductivity may be found easily. 
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C.2 Recursive Differentials 

Recalling from Appendix A that; 


■*'0 

R(>.) + i - r ■ 

K o + Z o 


(C . 9) 


and differentiating with respect to the layer thickness d^, 
we find 


Hf 1 - -t 0 <K 0 + V 


-2 


3? *0 

3d 1 


-a^K^ j (K 2 - Z 2 ) (sech 2 oudi) 


(K 0 +Z 0 ) 2 1 (Ki + Z j tanh uidi) 2 J 


(C.10) 


He can deduce immediately from equations (C.10) and (A. 4) 
that the general expression for 3Z^/3d^ +1 is of the form; 


3Z i 


5cT 


i+1 


“i+l K i+1 


& 


i+l’ - Z i + 1 !) <8eoh! a i+l d i+l 




<K i+l + 2 i+l tsnh “i+l d i+l> 2 


7 


(c.n; 

Noting that we can put the partial derivative of R(A) with 
respect to any layer (i+l) in the form: 

1 

3Z . 


3R( A) 

1 K 0 


1 

o 

N 


r-K 0 i 

id i + i 1 

L _ «W*J 

1 

L^J 


< W > J 


3Z 0 9Z X 


(C. 12) 

we see that a genera] form of the partials 3Z^/3Z.^ is 
necessary to allow a general expression for the partials 
of R( A) with respect to layer thicknesses. 
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Noting 
respect to 


that we can write the partial of 
in the form: 


Zq with 


9Z q K 1 z (1 - tanh 2 a ^d^) 
■ 

1 (K 1 + Z ^ tanh a^d^) 


(C . 13) 


and again inspecting equation (A. 4 ) , we can easily deduce 
the generalization of (C.13) as: 


9Z. 


K. 


i+1 


(1 - tanh 2 a. 


i+1 d i+i 


(K i+l + 2 i+1 tanh a i+l d i+l ) 2 


(C .14) 


Accordingly, we can write equation (C.12) in the general 
form: 


3R(1 ) 
3d i+l " 


- K, 


i-1 3Z . 

n 1 


3Z. 


(K 0 +Z 0 ) 2 ^0 3Z i+ ^ 


9d i+l 


- K, 


( w 2 


i K. 2 (l - tanh 

n — 


'°i a i> 1 


i=l (K^ + Z^ tanh a^d^) 2 ^J 


( a i+l K i+1 (K i+l 2 " 2 i+1 2 ^ sech 2q i+l d i +l I 
(K i+i + z i + i tanh a i+i a i + i> 


(C.15) 


> 


93 


ThuSf using the above recurrence relationships for the 
partials of R(X) with respect to d i? we can readily 
determine all the derivatives of d^. 

Recurrence relationships for the partial derivatives 
of R (X) with respect to the conductivities can be 
similarly obtained with the realization that is 
implicit in the characteristic impedances as well as 
the port impedances Z^, 

The solution posed herein should reduce the amount of 
computer time necessary for the inversion of data for the 
simple geologic case of one dimensional layering. However, 
more work must be done to determine the interdependence of 
0 . and d^ and the optimum iterative scheme for achieving 
data inversion. 
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Appendix D : The Image Approximation for a Multilayered 

Earth 

Wait (1969) and Weaver (1971) developed approximate 
solutions for the electromagnetic fields of a magnetic 
dipole over a conducting half-space. Earlier, Wait (1962 
and 1970) had developed expressions for a line source over 
a conducting layered medium. Additionally, Harrington 
'1973) suggested that the image approximation could be 
extended to the determination of the magnetic field of a 
magnetic dipole over a multilayered, conducting earth 
with the introduction of a "stratified image depth", alpha 
a. In this appendix, we shall outline the mathematical 
basis for Harrington's extension following the analyse:; of 
Wait (1962, 1969). We shall change our notation in this 
appendix to distinguish the ou of Chapter 4 from the image 
depth a. 

Consider a magnetic dipole of moment m at a height h 
above a multi-layered conducting earth spanned by the 
coordinate system (r,z) as illustrated in Figure D-l. 



Figure D-l. Magnetic Dipole Geometry. 


The quasistatic solution for the Hertz potential for a 
magnetic dipole can be written in the form: 


r = 


m I 1_ 

k \ R( 


+ P 


where 

P 


! u e " X . ,.-X(h 

j or +x J o (Xr)e 

o 6 


- z) 


dX 


R = (r 2 + (z - h) 2 ) 


1 

T 


(D.l) 


and X is a variable of integration and for M layers (Wait, 
1962) , u e is: 


u = u. 


Z 2^ K 1 + tan ^ u l^l 
1 + tan ^ u i^l 


where: 

u_ 


2 \ T 


m 


(X z + Y m ) 
m 


m = 0, 1 , 2 , 3* • • M 


2 

• = (ia u w - e u u 2 )= : i 0 u u) 

m nvm nr m 1 m M m 


for <^>>£^ 0 ) • 


(D.2) 


(D.3) 
(D . 4) 


96 




es 


TTfir-ini 



Vi + *m tanh VVi 


K + Z , tanh u h 
m m +1 mm 


(D.5) 


7 — v 

"M-l M-l 


K., + K.. , tanh u,, n h M , 

M M-l M~1 M-l 

K., i + k m tanh u„ h , 
M-l M M-l M-l 


K = 
m 


a + ioje 
m m 


U Y 

m _ j_m 

a o 

m m 


(D.6) 


and h is the thickness of the m — slab, 
m 

The major contribution to the integral P in (El) occurs 
for values of A .of the order of |y |. Thus, for |y m I > > I "Y Q I »• 
equation (D.3) becomes: 


u = (A 2 + y 2 ) 2 - 
m m 'm 


To this order of approximation 


(D.7) 


u e “ + ^e )2 a Y e 


(D. 8) 


(Z 2 - , ^ k i) + tanh u^h-^ 


Y e ~ Y i?l +(Z 2 /K 1 ) tanh u^^ 


(D.9) 


We now follow closely Wait 1 s (1969) analysis to determine 
the image approximation to the Hertz potential T . 
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Consider the Maclaurin series expansion of the 
expression : 


, aX u e ~ ^ “ .. 

fU> = e — p-y = £ V 

e m=u 

where a is to be determined and, 

a = ~ f (m) (0) 
m mi 

Carrying out the derivative operations, we find that 
. (u - X) 2 f 0 

f ' (*> - e — — y ~ ~ 

yi v e 


f' ' (A) = f' (A) \<x - + 2f (A)- 

t U eJ 

etc. By choosing a = 2/y g we find that a Q = 1 , 

a i =0 ' a 2 =0 ’ a 3^) 3 7i , a 4 =0 ' a 5 = " { 2 ) 18/5 ! ' a 6 = ^') 
etc. 

Thus: 


^74 - * 


(?) {01 + ••») 


and : 


P = + \ (1 + f fy 3 + 000 ) J 0 (^p)e 


-A (a - z + h) 




a 3 2 3 


V I OOO II 

™ ) S a 


where R = (r + (h + a - z) ) 
a 


(D, 10) 


(D . 11) 


(D, 12) 


(D , 13) 


40/6! , 


(D.ll) 


(D , 15) 


(D. 16) 


For c >>e a) and y A =u for all m, we can represent a as: 
mm urn 

/- h ^ 6 

iii ^ i 


tanh s(l+j)T” ) + Ri T~ 
a = 2y e = i_(l - j) j ‘-t- fc- 

1 + Rz tanh - (l+j)-j— ■ 


(D. 17 ) 


where: 


Ri 


t h a\ 5 1 

tanh u 1+ j ) & 2 j + R2^2 

" — ■ ——II — ■ 

1 + R2 tanh . ( 1+ j ) ^ 


2 J 


and: 




M+l 


^ ^ tanh ^ (1+j)^. + R m -j— 

M_1 1 + RjJ *Wl tanh ~ ( 1+ j ) ^ 


M 


J M 


(D. 18 ) 


(D . 191 


For the two layer case, h 2 -*■ • and 63 6 2 . Thu. p ; 


Ri 


1 + R 2 = , 
1 + R z 


(D. 20 ) 


and : 


a = 6 1 (1-j) 


tanh j<l+j> + |f 

1 + §7 tanh *' {1+j) IT t 
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As shown in Chapter 4, we can write; 


I 

! 

I 

I 

i 

I 

I 

I 

I 

l 


H = 


a z r 

3z r 


and 



(D.21) 


(D. 22) 


Retaining only the first term of the Maclaurin series 
expansion in equation (DlIS) and normalizing the magnetic 
fields by multiplying H and H by 4irr 3 /m, we find 
expressions for the normalized "image" fields near the 
surface to be: 

H ~ 3z 3 (a/r) + H.O.T. (neglected) (D.23) 

r * r U + (a/r) 2 j s / 2 

H s 1 - + H.O.T. (neglected) (D.24) 

2 U + (a/r) 5/2 

Noting that the higher order neglected terms contain 
expressions of the form (I) : 

00 

I = 2 8 (1 + <a/r) 2 )” m/2 (D.25) 

n>5 m 


we can readily deduce that equations (D.23) and (D.24) are 
accurate descriptions ot the magnetic field intensities 
only for ranges large with respect to the real part of a, 
i.e. r r>>Real a. Graphically comparing the expressions -or 
the "image" fields with those of the exact integral fields 
for the H component (Figure d- 2) , we note that, when the 
source is on or near the surface of the earth, the image 
formulation for the magnetic field near the surface is 
accurate within 10% only for r > 7 Re|a| and is not of much 
practical use in our application. 
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Figure D-2. Comparison of H for Integral 
and Image formulations. 
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Appendix E: Switched Source Synthesizer 

A practical implementation of the SWEEP inductive 

coupling geophysical prospecting technique requires that 

the source magnetic field be generated with a physically 

small coil that can be deployed easily or, preferably, can 

be mounted permanently on a small vehicle. But, to permit 

examination of strata to depths of the order of a few 

hundred meters we must use a source with a magnetic moment 

4 2 

of the order of 10 amperes-turns-meters . Practically,, we 
need currents of the order of 50 amperes and power of the 
order of 15 kilowatts. 

For some time now work has been proceeding on the 
development of an all solid state high power DC-to-AC 
inverter for another application which also happens to 
require currents and powers of this order. And this work 
has led to the development of novel techniques^' which 
allow hnhly efficient syntheses of various waveforms, 
including sinusoids, directly from storage batteries. 

The concept for the switched source synthesizer is 
shown in Figure E-l. The configuration allows various 
combinations of batteries to be connected sequentially in 
series to generate a time varying voltage. In the 
illustration, the battery voltages are V^/2, V /3 , and 
V /6, respectively, and when the batteries are connected 
in a proper time sequence, the arrangement 
generates a quasisinusoidal waveform as shown in 
Figure E-2. 


(1) R.H. Baker and L.H. Bannister, "Electric Power Converter", 
Patent Pending. 
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Figure E-3. 


Harmonic Content of 
Six-*tep Sinewave 
Approximation 
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There are several interesting technical features of a 
switched source synthesizer that make it quite attractive 
as a DC to AC inverter: 

a. For the sinusoidal waveform shown in Figure E-2, 
not only is the total harmonic distortion small 
(6%) , but in addition, the major components of 
distortion are concentrated at the higher harmonic 
frequencies. Figure E-3 shows the relative 
magnitude of the various harmonic components for 

a six -step sinusoidal approximation. The largest 
component of distortion (3%) occurs at the 35^ 
harmonic which is easy to filter from the 
fundamental. 

b. Any frequency from near zero to several 
kilohertz can be generated with constant 
amplitude and harmonic distortion because the 
inverter does not use a transformer. 

c. The switched source synthesizer is inherently 
efficient because it is a nonlinear system and 
the transistor switches are never operated in 
the linear region of simultaneous high voltage 
and high current. This high power efficiency 
(n>95%) allows circuit packaging in compact 
form which in turn leads to light weight 
portable equipment. 

We believe that the switched source synthesizer 
configuration combines the features of portability, high 
efficiency, and low harmonic distortion over a wide 
frequency range, which make it an almost ideal source for 
our SWEEP inductive coupling geophysical probing technique. 
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